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Abstract Measurements of the order parameter texture of rotating superfiuid 
3 He-B have been performed as a function of the applied azimuthal counterflow 
velocity down to temperatures of 0.2 T c . The results are compared to the hydro- 
static theory of 3 He-B. Good agreement is found at all measured temperatures 
and rotation velocities when the flow anisotropy contribution to the textural free 
energy is adjusted. This gives a superfiuid energy gap A(T) which agrees with 
that measured by Todoshchenko et al, with A(0) = 1.97fc B T c at 29.0 bar. The 
B-phase susceptibility, longitudinal resonance frequency, and textural phase tran- 
sition have been extracted from the measurements as a function of temperature 
and azimuthal counterflow velocity. Owing to decreasing absorption intensities the 
present measuring method, based on the line shape analysis of the NMR spectrum, 
loses its sensitivity with decreasing temperature. However, we find that in prac- 
tice the measurement of vortex numbers and counterflow velocities is still feasible 
down to 0.2 T c 

Keywords Energy gap • Density anisotropy ■ Flare out texture • Nuclear 
magnetic resonance ■ Superfiuid counterflow 
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1 Introduction 

Non-invasive nuclear magnetic resonance (NMR) measurements on superfiuid 3 He- 
B in a rotating cylindrical environment provide information on the axially sym- 
metric order parameter texture. This method - based on an analysis of the line 
shape of the NMR spectrum - has been the foremost means of measuring the 
number of rectilinear vortices in a stationary state of rotation [I] . 

A prominent feature in the NMR line shape is the so-called counterflow (hence- 
forth cf ) peak, which appears when there is sufficient flow velocity in the azimuthal 
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direction. Counterflow is the difference between the velocities v n of the normal 
component and that of the superfluid component v s in the presence of A" axially 
aligned vortices, i.e. v c f = v n — v s . In stationary state, rectilinear vortices form 
a cluster along the axis of rotation with a distribution which corresponds to the 
equilibrium vortex number density of the applied rotation. The radius of the clus- 
ter can be smaller than that of the system if A" < N v , where N v w 2itR 2 f2 v j 'k is 
the equilibrium number of vortices. Here n is the quantum of circulation, R the 
radius of the system and f2 v the rotation velocity associated with A" v vortices in 
equilibrium state. For an example of a vortex cluster configuration, see Fig. 1 in 
[2]. Inside the vortex cluster the average superfluid velocity (v s ) w v n = Qr, while 
outside the centrally located cluster u s = nN/2irr. The local cf velocity v c f is thus 
approximately zero inside the cluster, while outside 



In contrast, flow in the axial direction is expressed by an increase of absorption 
in the spectrum near the Larmor peak. With continuous wave (cw) NMR, the 
sweep time of the whole spectrum can be significantly longer than the characteristic 
time of the dynamics, in which case one has to resort to sweeping that part of the 
spectrum which expresses the vortex dynamics best: absorption of the cf-peak or 
absorption near the Larmor frequency. In this paper we only consider states where 
transient effects have ceased and the vortex configuration with rectilinear vortices 
is in the stationary state. 

Examples of NMR spectra, measured in different conditions of flow at 0.2 
0.3 T c temperatures, are shown in Fig. [1] The spectra display absorption peaks 
which are shifted differently from the Larmor value. The shift and peak intensity of 
these absorption maxima identify the flow conditions. The integrated absorption 
intensity and simultaneously the shift of the peak at given rotation velocity (2 
specify the number of vortices N in the central vortex cluster and thereby the 
azimuthal counterflow circulating around the cluster. 

This NMR measurement technique based on the line shape analysis has been 
extensively used, e.g. in measurements on dynamic remanent vortices [3J, the onset 
temperature T on to turbulence [3], turbulent vortex front propagation [5] and the 
stability of flow in the low temperature regime [6 . While in some experiments the 
measured spectra are used in a qualitative way (e.g. the presence or absence of the 
cf-peak), measurements where the number of rectilinear vortices in the sample is 
important need an (experimental) calibration for the number of vortices (or the 
azimuthal cf expressed as rotation velocity Q c { = fl — J? v ) as a function of the cf 
peak height. 

In this paper we compare the measured NMR response with the calculated 
spectra quantitatively. We start by reminding about the hydrostatic model of 
superfluid 3 He-B in a cylindrical geometry 0(8] • We discuss the superfluid order 
parameter textures, textural energies in applied flow and the relation to NMR. Af- 
ter the discussion of the experimental setup, we show results on measured texture 
transitions, on the susceptibility, and list spectra as a function of the azimuthal cf 
velocity. We compare measurements with calculations and obtain good agreement 
when we fit the density anisotropy, Ajjv, of the superfluid. The discussion on the 
measured parameter Ahv includes its relation to the energy gap A and consid- 
erations about its magnitude. We conclude the paper with the identification of a 
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Fig. 1 Order parameter textures. The measured and calculated NMR line shapes of 
the simple, parted and extended flare out textures of 3 He-B in a cylindrical cell with radius 
R = 3 mm. The measured spectra are presented as a function (horizontal axis) of the magnetic 
field H, or equivalently the frequency shift Av = v — v\ J from the Larmor frequency ^l, or as a 
reduced frequency shift Av = sin 2 /3. Here the line shapes are normalized to the latter axis. The 
associated textures, i.e. n- vector representation (by means of the angles a and (3) as a function 
of reduced radius are depicted in the inset and on the left. Parameters of the simple flare out 
texture: fork temperature T = 0.25 T c , pressure p = 29 bar, rotation velocity Q = 0.2rad/s 
with no vortices, cw-NMR frequency v r f = 1.967 MHz and Larmor field = 60.65 mT; 
parted flare out texture: T = 0.31 T c , f2 = 0.9rad/s with vortex cluster S7 V = O.lrad/s, 
u Tl = 965.0kHz and _ff L = 29.75 mT; extended flare out texture: T = 0.36T C , Q = 0.7rad/s 
with no vortices, v T f = 965.0 kHz and = 29.75 mT. In the calculations, temperature T 
(controlling the temperature dependent frequency shift Qb) was a fitting parameter. Other 
fitting parameters; simple flare out texture: field velocity parameter Ajjv = 1.4kg/m 3 T 2 and 
field inhomogeneity AH/H = 8.3 ■ 10 — 4 ; parted flare out texture: Ahv = 3.1kg/m 3 T 2 and 
AH/H = 8.8 ■ 10~ 4 ; extended flare out texture: A HV = 5.9kg/m 3 T 2 and AH/H = 8.8 ■ 10~ 4 . 
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newly found resonance peak in the NMR spectrum. The peak is formed due to a 
non-local resonance in the potential well created by the texture near the cylindrical 
boundary of the experimental cell. 



2 Hydrostatic Theory of 3 He-B 

In the B-phas^3 of 3 He, the Cooper pair with total spin s and angular momentum 
1 is condensed to a state where the angle 9 between s and 1 is fixed, due to 
interactions between the nuclear dipole moments of the helium atoms, to 9 = 
arccos(— 4) w 104° [9]. In this configuration the total angular momentum j = 0, 
while 1 = 8 = 1 [10] . The order parameter matrix can be written as [11] 

A^ = A„ v R vj (h,9)e^, (2) 

with A(T) the temperature dependent energy-gap, </> the condensate's phase, and 
Rij the rotation matrix describing the rotation of the spin and orbital coordinates 
relative to each other by 9 around the axis oriented along the unit vector n: 

Rij(n,9) = cos 9 5ij + (1 — cos9)fiihj — sm9 eyfcfife. (3) 

Spherical coordinates are used for the unit vector n with a and (3 denoting the 
azimuthal and polar angles with respect to the axis of the static magnetic field H. 
The orientational distribution of the n- vector over real space is called the texture. 
On length scales longer than the dipole-dipole interaction £d ss 10 /im, a number 
of different interactions give rise to textural anisotropy [13_ in the bulk and at the 
surface. We start with four bulk terms and continue later with the surface term. 
The magnetic orientational free energy term 



-a J d 3 r(n • H) 2 (4) 



results from the gap anisotropy in a magnetic field H . The dipole velocity term is 
Fov = -A DV f d 3 r[n • (v. - v n )] 2 . (5) 



The gradient term 



Fq = I d 3 r 



dRai dRaj , ORaj ORaj 
A G1— ^ a r A G2 — ^ 7^ 



(G) 



reflects the coherence of the superfluid state - the stiffness of the order parameter 
with respect to spatial change - and resists rapid spatial change in the texture. 
The field velocity term 



-At 



- 1 d 3 r[U ■ It ■ (v. - v n )] 2 (7) 



arises from the anisotropy axis 1 preferring to orient in the direction of the flow, 
which changes the direction of the anisotropy axis of the energy gap A. 



2 The B-phase is also referred to as the BW-state after the discoverers R. Balian and N.R. 
Werthamer [10] 
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We continue with the surface term. Here we assume that the curvature of the 
surface boundary is small and that the length scale of the distorted region is small 
compared to the dipole length £d . The surface gives rise to the surface field term 



F S H 



-d 



[ d 2 r(H 
Js 



R 



(8) 



Here the unit vector s is perpendicular to the surface and points toward the su- 
perfluid. For explicit expressions of the parameters a, Adv, Aqi, ^G2, d we refer 
to formulas (35)... (42) in [7]. The field velocity parameter Ahv is discussed here 
explicitly, since we use it later on as a fitting parameter, and to extract information 
on the energy gap A(T): 



<HV 



a 2 (i + ±Ffy) 2 
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(9) 



with p the fluid density, m* the effective mass, ^/2tt = — 32.435MHz/T the gyro- 
magnetic ratio, Ff the symmetric and F ( a the anti-symmetric Fermi-liquid param- 
eter. The temperature-dependent functions Zj are defined as 



-J-/2 



(10) 



where the Matsubara GIlGrglGS cirG En = nT(2n - 1) with n = 0, ±1, ±oo and 
the Yosida function is given by Y = 1 — Z 3 (T). 

So the question arises: how do the bulk terms (HJ)...© and the surface term 
(|8| affect the texture in a cylindrical environment? 

The rotational symmetry and continuity at r < R means that at r = 0: n || H, 
i.e. p = 0. At the boundary, the orbital momentum density L is perpendicular to 
the surface and sin 2 f3 = 0.8 (i.e. j3 « 63.4°). In the intermediate region the texture 
is smoothed out by the gradient term Fq. When the superfluid is stationary, i.e. 
non-rotating, the simple flare out texture is the lowest free energy state. See Fig. 
[T]for an example of the simple flare out texture. 

By rotating and applying azimuthal counterflow to the system, a negative 
energy contribution -Fhv arises from the 1 vector orienting along the direction of 
the flow, which changes the direction of the density anisotropy, Sp^ = pjj"^ — Pj" , 
and reduces the total energy. This has a direct effect on the /3 angle, see inset in 
Fig. [T]of the parted flare out texture. The ends of the /3 curve stay fixed at r = 
and r = R, but the intermediate slope moves to the left. At the critical velocity 
Q c \ the curve starts to show a plane near the surface boundary (where the cf is 
the largest). The angle f3 on the surface is slightly moved due to the stress from 
the surface gradient energy term Fsg . The resulting texture is now the parted flare 
out texture with the cf-peak as a growing prominent feature. 

At the critical velocity f2 C 2, the region with almost constant (3 has grown to 
its maximum, but not minimal in energy due to the surface gradient energy term 
Fsg- The texture can lower its energy by flipping the angle f3 at the surface and 
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the transition to the extended flare out texture occur sE Here /3 now passes across 
90° in the vicinity of the surface boundary. Further increase of the applied flow 
extends the region with a near constant j3 even more, until saturation. 



3 Nuclear Magnetic Resonance 

The details of the texture would be hidden from observation were it not for the 
fact that spin-orbit interaction with the magnetic moment of the 3 He nucleus allow 
us to probe the texture non-invasively pQ . The Leggett equations 

86' 

?L = - 7 H to t + ^-s (11) 

Ot XB 

5 = 7 S x H tot - ^Mn(n • 0') (12) 
ot 7 2 

describe the response to small rotations 6' of the spin density S, where Xb(T) is the 
temperature dependent B-phase susceptibility, fis (T) the temperature dependent 
B-phase longitudinal resonance frequency and Htot the sum of the static field 
and a small transverse rf field oc e lujt . In the high- field limit (H 3> 2.5 mT), the 
transverse resonance frequency reduces to 

p 2 

lo'I + f2 B sin 2 P « o; L + sin 2 /3 (13) 
2wl 

with cjl = |7|-ff the Larmor frequency. Here we primarily look at the frequency 
shift from the Larmor value. In continuous wave (cw) NMR experiments, the 
NMR response is measured as a function of the longitudinal magnetic field H. 
The spectrum is then converted to the frequency domain using the relation 

Av= iJh ± _H)H VrU (14) 

where u r f is the frequency of the transverse excitation field, H is the measured 
field and Hi, the Larmor field. The latter can be obtained from the NMR response 
in the normal phase. We define the reduced frequency shift Av using (1131) as 



2w L , 



Av = 7^(^ - wl) = sin 2 p. (15) 

The Kramers-Kronig relation tells that the absorption and dispersion signals 
are interrelated and the area below the measured absorption line shape xty) is 
a constant. Measured NMR line shapes can thus consistently be compared with 
e.g. calculations in different domains: in the frequency domain v or in the reduced 
frequency domain v. The normalization of the spectra is done by dividing with 
the integrated area under the line shape above the baseline: 

m = t4^p- d7) 



1 There is ambiguity in the terminology in different publications: a) flare out and flare in 
is sometimes used instead of parted and extended flare out; b) flare out sometimes refers to 
both parted and extended flare out texture |3l ll4|[2l) 
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The normalized NMR line shape x m the reduced frequency domain Av is free of 
the NMR setup parameter v r f . For examples of conversions between the magnetic 
field H, frequency shift Av and reduced frequency shift AD, see the different 
horizontal axes of the NMR line shapes in Fig. [1] 



4 Numerical Calculations 

Numerical calculations have been performed by finding the equilibrium texture, 
i.e. the texture which corresponds to the minimum configuration of the free energy 
and by determining the corresponding distribution of NMR intensity [8] . The cal- 
culations use the following input: for geometrical parameters the cylinder radius 
R, rotation velocity fl and the number of vortices N, which for convenience we 
characterize by the rotation velocity i? v , at which this number of vortices is the 
equilibrium amount; for the environmental parameters the reduced temperature 
T/T c and pressure p is used; for the NMR parameters the Larmor frequency 
and the inhomogeneity of the applied magnetic field AH/H; and for the physical 
parameter the field velocity parameter Ahv- Other physical parameters are cal- 
culated directly from theory using temperature T and pressure p, or taken from 
experiment (see Fig. 1...4 of [7]). 

The outline of the calculation is as follows: from the geometrical parameters the 
vortex configuration is determined and the counterflow profile v c f (r) is calculated. 
Using ([I]) . . . (0) together with the environmental parameters, the textural energy is 
minimized using the truncated Newton algorithm 12j which gives the solution as 
the n-vector expressed in a(r) and j3{r). The discretization step for r is 15/im. 

For the determination of the NMR absorption spectrum, the so-called local 
oscillator model is used: the fluid as an assembly of uncoupled oscillators with 
frequencies determined by the local value of /3(r), see (|13p . The NMR line shape 
is then a sum of the individual contributions in the volume V 



The local inhomogeneity of the magnetic field AH/H causes dephasing on the 
local oscillators due to the spread AH around the average field value H. The 
characteristic time th is of order 



By introducing the term — ThS^ on the right-hand side of (|12p. where Sj_ is the 
transverse component of the magnetization, the delta function in (|18|) becomes 
effectively a Lorentzian with a width Alu = Th- We determined the average field 
inhomogeneity for Aui by measuring the NMR response in the normal phase. The 
superfluid texture is of the flare out type: the angle (3 changes as a function of 
radius r, hence the measured NMR line shape of the flare-out texture depends on 
the field inhomogeneity at r. In combination with that the NMR pickup coil does 
not probe the fluid homogeneously, it makes the field inhomogeneity as measured 
in the normal phase only indicative. Therefore, the field inhomogeneity parameter 
is a fitting parameter when measured and calculated spectra are compared. 




(18) 




(19) 
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u |^ 6 mm 




-1 o 1 with orifice 



r (cm) 

Fig. 2 Experimental setup in 2005: the left figure shows the experimental quartz cylinder 
with a detector coil at both ends of the NMR sample volume. The bottom spectrometer was 
used for the linear NMR (full spectra) measurements, while the top spectrometer was used for 
non-linear NMR (magnons) measurements. The field inhomogeneity was determined by fitting 
the superfluid NMR line shape with calculations and was consistently AH/H si 8.8 ■ 10 — 4 at 
all measured temperatures. The right figure shows schematically the LC-resonance circuit and 
the NMR Q-meter setup. The superfluid is excited with a high-Q LC-resonator; the magnetic 
pickup changes the inductance which is measured by the pickup loop. The signal is compensated 
by a second phase locked function generator for high resolution measurement after amplification 
using the Stanford Research Systems SR560 differential amplifier. A phase lock-in amplifier 
(Stanford Research Systems model SR844) is used for data acquisition. 



The first order texture transition from the parted to the extended flare out 
texture and back, shows hysteresis in the dependence on f2 and T, therefore dual 
solutions exist for the texture around the critical velocity Q C 2- To cope with this 
degeneracy, the program accepts a(r) and /3(r) as initial starting values and then 
proceeds to find the closest local minimum of the textural free energy. The initial 
guess can be as crude as a straight line from /3(0) = to /3(R) = 63.4° or /3(R) = 
180° — 63.4°, for the parted and extended flare out textures respectively. 



5 Experimental Techniques 

NMR measurements were performed with two different setups, which both in- 
cluded a long cylindrical container with radius R = 3 mm mounted on the nu- 
clear stage of a rotating cryostat. The cryostat rotates up to 3rad/s with an 
^-dependent noise level which is Afl = 0.006 rad/s (peak-to-peak) at lrad/s. 

The 2005 setup (referring to the starting year) had above the sinter on top of 
the nuclear stage two quartz tuning fork oscillators (henceforth forks). The main 
sample volume is separated from the fork volume by a plate with an orifice (0 = 
0.75 mm). The plate prevents vortices from entering the main volume. At each end 
of the main volume, two NMR pick up coils were mounted for continuous wave 
(cw) NMR. See Fig.[2]for the cell and the NMR resonance circuit. For the flare out 
texture measurements in this paper, we use the bottom spectrometer. The order 
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parameter texture was probed using a resonance circuit tuned to v T f = 965.0 kHz 
with a quality factor Q = 6050 and a Larmor field of Hi = 29.75 mT. 

The 2009 setup had two forks in an isolated volume, which sits between the 
sinter volume and the main volume. Two orifices connect the three volumes. The 
orifice on the main volume side was bigger than the orifice on the sinter side, 
such that at low temperatures the quasi-particle flux from the main volume side 
dominates the quasiparticle density of the fork volume [16] . The NMR coil at the 
bottom end was used for full spectra measurements. The resonance circuit had 
a pre-amp at 4K and a second amplifier at room temperature. The resonance 
frequency of v Y i = 1.967 MHz was used, which corresponds to the Larmor field 
Hi = 60.65 mT. The resonance circuit had a quality factor of Q = 3900. The 
transverse rf field of the spectrometer coil was H r f = 3.4 nT and the coil had an 
inductance of L = 12.7/^H. 

In both setups the temperature was measured using a fork, which was cali- 
brated above 0.35 T c against a Melting Curve Thermometer (MCT) mounted on 
the nuclear stage [TT] . In the low temperature regime where the MCT loses ther- 
mal contact with the 3 He NMR sample, the resonance width of the fork response 
was interpreted to be proportional to the quasi-particle density which depends 
exponentially on temperature [11] . The heatleak was determined to be Q = 20pW 
(at J? = lrad/s) for the 2009 setup, while for the 2005 setup this was estimated 
to be of one order of magnitude larger. 

The superfluid was continuously probed by sweeping slowly the magnetic field 
in the NMR magnet and measuring the absorption x(^) an d dispersion in the 
NMR pickup coils at 4 samples/s with a lock-in amplifier. The typical sweep rate 
of the main NMR field was H = 9 -r 13/iT/s. 

The magnetic field of the NMR is produced by a combination of solenoid coils 
which by design have an axial and radial inhomogeneity of AH Z /H Z = 7.3 ■ 10 -5 
and AHr/Hr = 2.0 • 10 -5 , respectively. Niobium shields protect against external 
magnetic influences. In practice the field inhomogeneity is increased by an order 
of magnitude by the presence of the superconducting NMR pickup coils. The field 
inhomogeneity was measured a number of times (after different cool downs to liquid 
helium temperatures) in the Fermi-liquid state (normal phase) with low enough 
excitation not to saturate the signal. The NMR response was distorted and far from 
an ideal Lorentzian. A poor fit gave an inhomogeneity of 1.7 • 10 -3 . Measurements 
of the field inhomogeneity were taken with different magnetic fields -ffdemag = 
0.1 -r 1.1T in the superconducting magnet used for adiabatic demagnetization 
cooling (or nuclear cooling), to verify that the superconducting Nb shields around 
the NMR setup are protecting the NMR measurement from external magnetic 
influences. The discussion on the field inhomogeneity is continued in Section [8] 
where we use the inhomogeneity as a fitting parameter. 

The NMR line shapes were processed by substracting the baseline and correct- 
ing the phase using both the measured absorption and dispersion signals of the 
lock-in amplifier. The magnetic field sweep was linearly corrected for the inductive 
delay of the magnet at sweep rate H\ this correction made up and down sweeps 
to overlap. 

Since we want to measure the NMR line shape of the vortex free state as well as 
those with different number of rectilinear vortices N, the sample state is prepared 
above the onset temperature T on of turbulence [3l|3j . A rotating vortex free state 
is created by increasing the rotation velocity fl starting from zero with Q = 
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0.03rad/s 2 , at a temperature T = 0.7-^0.75 T c . This temperature is above the onset 
temperature T on to turbulence and below the AB-phase transition temperature 
Tab- To prepare the superfluid with a predefined vortex cluster, the following 
recipe was used: at these high temperatures the superfluid is rotated above the 
critical rotation velocity J? cr (in our cell fi C i = 2.5 -f- 3.5rad/s) where vortices are 
created when the cf velocity reaches the bulk liquid critical value at the largest 
imperfections on the cylindrical surface. At these high velocities the cf-peak drops 
significantly indicating that a large number of vortices are created. The rotation 
is then slowed down to the desired rotation velocity fi v . After annihilation of 
extra vortices and reaching the equilibrium vortex state, the cylinder contains 
iV v rectilinear vortices. The rotation velocity is then increased to a desired value 
fi < fl cx making sure that no new vortices are created. The cell is now rotating at 
a rotation velocity fi with a known number of vortices N, which is smaller than 
the equilibrium number of vortices at this rotation velocity [6 . 

The measurements in the 2005 setup were performed by recording the NMR 
line shape during a slow cool down (T « 50 • 10 T c /s) to low temperatures. The 
cool down was frequently interrupted for temperature stabilization/verification. 
The rotation velocities used were fi = 0.7, 0.8, 0.9, 1.05 rad/s, with clusters of 
vortices corresponding to J? v = 0, 0.1, 0.25rad/s. These low numbers of vortices 
have a significant effect on the texture pQ. After reaching low temperatures, the 
temperature was swept back to the starting temperature to compare the initial 
and final NMR line shapes to make sure that no new vortices had been created 
in this temperature cycle. Note that this cycle takes place to a large extent in the 
regime where superfluid 3 He-B can have turbulent dynamics: below T ~ 0.6 T c 
the superfluid Reynolds number Re a increases rapidly from unity to Re a = 10 3 
at T = 0.2 T c |18| . However, the experiments presented in this paper are not 
affected by turbulence, due to the fact that there are no vortices connected to the 
cylindrical boundary [6] . 

In the measurements taken in the 2009 setup, the cell was prepared in the 
vortex free state by increasing the rotation velocity at T w 0.7T C without the 
creation of vortices. The cell was then cooled down at constant rotation fii to a 
preselected temperature T where the NMR line shapes were measured in the vortex 
free state at different rotation velocities fi in the range 0.35 H- 2 rad/s. Between 
the measurements at the varying rotation velocities fi, the rotation was brought 
back to the initial rotation velocity fii for verification that no new vortices had 
been created by comparing the initial and current NMR line shapes. Temperature 
control was stable within AT = 0.005 T c , as measured by the forks. 



6 Dynamic Susceptibility Measurement 

We use the area of our spectra as a measure for the static susceptibility by com- 
paring the absorption x( Cl ") m the superfluid state with the absorption xn in the 
normal state above T c to assure linear NMR was measured. The signal amplitude 
relates to the susceptibility xb{T) as 



xM _ ^oTtXbCT) Vs{uj) 
Xn 2 xn / V B (w')dw' ' 



(20) 
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Fig. 3 Magnetic susceptibility of 3 He-B: The reduced magnetic susceptibility Xb/xn as 
a function of reduced temperature T/T c . The results are obtained by taking the ratio between 
the integrated (baseline, phase and amplitude corrected) NMR line shape in the vortex free 
supcrfluid state (Landau state) and the Fermi-liquid state (normal state). Measurements are 
shown taken at different Q (gray markers) in the range 17 = 0. 15-^2. Olrad/s. The scatter results 
from the baseline drift during the measurement. There was no cross-correlation between the 
scatter in Xb/xn and Q, i.e. random scatter around the average. Averages (blue markers) are 
taken on isothermals with the error-bars indicating one standard deviation from the average 
value. Our measurements are compared with measurements by Webb et al. 121 1 and by Scholz 
1221 published in reference 1231 . The two curves represent the theoretical model fl22[) with 
different values for the energy gap A which enters through the Yosida function. The solid line 
shows the model with A BC (T) and the dashed line with a scaled energy gap ^scaled CO ( see 
lt24l) and lt25l) for definitions). 



where Xb (T) / Xn is defined as the ratio of the total integrated NMR absorptions in 
the superfiuid phase and the normal phase which can be experimentally determined 

Xb(T) _ UV a (u/)dw') T f2n 

X N (JV s (0j')dcj')T o ' ( ' 

The results of xb{T) measurements of the vortex free state at different rotation 
velocities 1? as a function of temperature T are plotted in Fig. [3] The values of 
Xb at different fi have a larger scatter, but there was no correlation with the ap- 
plied velocity Q. During the measurements the baseline was drifting, which affects 
the normalization and is seen as the main cause of the scatter. The theoretical 
expression for the susceptibility according to the weak coupling theory is [20 

XN " U + t0 > 1 + F S (I + iy(T)) - 3 + 2F§ ■ (22) 
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At a pressure of 29bar, -F a (29) = -0.75, which gives xb/xn = 0.33. 

In Fig.|3] xb/xn is plotted for the energy gaps A SC (T) and Z\ scaled (T). The gap 
Z\ sc is the weak-coupling gap Z\ WC (T) with corrections using the strong-coupling 
approximation. The scaled gap A sca ,i e d{T) is the strong- coupling corrected gap 
A sc linearly scaled such that it coincides in the low temperature limit with the 
experimental observed gap A exp as measured by Todoshchenko et al. [19]. To 
summarize: 

AvjciT) = gap weak — coupling, (23) 

^sc(^) = A WC (T^ (strong — coupling corrected, (^^) 

A scaled (T) = A sc (T)^^-. (25) 

The low temperature limiting values of the gap at p = 29 bar are Zi wc (0) = 
1.79fc B T c , Z\ sc (0) = 1.87fc B T c [7 and Z\ eX p(0) = 1.97fc B T c . 

The plot includes high temperature data from Webb et al. [21] and Scholz [22] 
at a variety of pressures. Both from their data and ours it is seen that the pressure 
dependence of Xb/xn, which comes about via Fq and the strong-coupling correc- 
tions, is not very strong above p = 20 bar. Also included is data from Ahonen et al. 
[28] which is taken at the same pressure (29 bar) as our measurement. The determi- 
nation of the static susceptibility of superfluid 3 He using the imaginary component 
of the dynamic susceptibility measurements (NMR) integrated over all frequencies, 
has earlier been in conflict with "true" static SQUID-based susceptibility measure- 
ments. Hahn et al. [23 suggest the difference with the NMR measurements, where 
the integrated area is used to determine the susceptibility, the integration range 
is limited to frequencies around Larmor i^l and the SQUID measurements include 
an extra contribution caused by, e.g. the electrons. For now, our measurements 
(at T < 0.7T C ) of the static susceptibility using NMR is consistent with the NMR 
data of Ahonen et al. and Scholz. 



7 Textural Transitions and CF-Peak Behavior 

In measurements at fixed temperature, we have observed textural transitions by 
changing the rotation velocity fl by a step (rapid acceleration at fl = 0.030 rad/s). 
The critical velocity fi c i of the second-order transition from the simple to the 
parted flare out texture is difficult to determine exactly since the transition is 
continuous. 

Fig. |31 shows the cf-peak height x' c { P (plotted in the reduced frequency domain) 
as a function of temperature for the parted flare out texture with no vortices 
present at several rotation velocities fl. At \ ~ 1-5 is the transition line where 
the cf-peak vanishes. What is left is the absorption of the simple flare out. The 
transition fl c \ does not seem to have any hysteresis within the resolution of the 
step Afl = ±0.050rad/s we used. 

In contrast, the second texture transition at fl C 2 between the parted and 
extended flare out textures has a large hysteresis in the i?-dependence and is 
of first-order. This transition happens instantaneously within our sampling rate 
of four samples/s. The hysteresis increases with decreasing temperature: around 
T = 0.5 T c the hysteresis is within Afl = 0.3 rad/s, while at T = 0.25 T c the 
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Fig. 4 Counterflow peak behavior near Q C V cf-peak height x'rfp plotted as a function 
of reduced temperature given in terms of the rotation velocity 17 in the vortex-free state. In 
the plot the cf-peak starts to appear at x' c { p ~ l-^i where the texture transforms from simple 
to the parted flare out texture. Lines connect data points with the same applied flow. With 
decreasing temperature, the measured cf-peak height decreases rapidly. 



hysteresis can exceed Afi = 0.8rad/s. Typically, the texture change seems to hap- 
pen during the step change in rotation velocity and rarely occurred during stable 
rotation. 

Fig.[5]shows the hysteresis of L? C 2 in the velocity dependence at 0.25 T c with a 
vortex cluster of L? v = 0.12 rad/s. A transition from the extended to the parted flare 
out texture occurs around J? « 0.75 rad/s when the velocity is step wise decreased. 
If the velocity is increased the texture change is reversed to the extended flare out 
texture at Q = 1.35 rad/s. 

In the measurements at constant L? the temperature was swept. The texture 
transition at J? C 2 between the parted and extended flare out textures is also ob- 
served to be hysteretic with a transition regime of AT = 0.025T C width. Fig. [6] 
shows the cf-peak behavior in a temperature cycle from high to low temperatures 
and back. In the low temperature regime, the temperature at which the transition 
from the extended to parted flare out texture occurs on a cool down is lower than 
the temperature of the transition from the parted to the extended flare out texture 
during a warm up. The large increase in absorption of the cf-peak from the ex- 
tended to parted flare out texture during cooling is due to the j3 angle not crossing 
the 90° angle in the parted flare out texture. The smaller /3 angle enhances the 
absorption of the cf-peak since the absolute value of the free energy term |.Fhv| 
is now smaller, and less energy is required to orient 1 along the direction of the 
flow. In effect 1 orientates along the flow over a larger region of space. Or in dif- 
ferent words: the cf peak increases in the texture change from the extended to 
parted flare out texture, since the spectrum width becomes smaller. The argument 
is reversed for the change in cf-peak absorption in the texture transition from the 
parted to the extended flare out texture: the angle /3 now crosses 90° and is on 
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Fig. 5 Velocity hysteresis of the critical texture transition fi c 2 between the parted 
and extended flare out texture: the plot shows the cf-peak at T = 0.25 T c with a vortex 
cluster of fi v = 0.12 rad/s during a velocity sweep down starting at 1 rad/s and up. In this 
cycle, the NMR spectra were recorded and the cf-peak was extracted after each spectrum was 
normalized in the reduced frequency domain. A transition occurred from the extended to the 
parted flare out texture when the velocity fi was lowered to fi = 0.8 rad/s. When the velocity 
was lowered further, the texture did not change. After step wise increasing fi, the texture 
stayed in the parted flare out texture until fi was increased to 1.35 rad/s where the texture 
flipped back to the extended texture. 



average steeper over the same space and the applied flow affects the tilt of the /3 
angle less. 

In Fig. [7] the textural transition Q c \ between the simple and parted flare out 
textures, as well as the textural transition J? C 2 between the parted and extended 
flare out textures are plotted as a function of reduced temperature. At high tem- 
peratures down to T w 0.4 T c , the transition velocity fl c \ is approximately con- 
stant. At temperatures below T w 0.4 T c , fi c i increases rapidly. The hysteresis in 
the transition velocity fi C 2 is approximately constant (AQ = 0.1 rad/s) down to 
« 0.4 T c where it starts to increase rapidly with decreasing temperature. Below 
0.25 T c , the hysteresis in applied counterflow is so big that once the system is in 
the parted flare out texture a transition back to the extended flare out texture was 
not observed below a critical velocity fi cr « 2.1 rad/s at which vortices nucleated 
due to imperfections of the cylindrical surface or multiplication of trapped vortices 
in the corner of the cylindrical container. Our result is compared with data ob- 
tained by Korhonen et al. [23] at similar magnetic fields H , but at differing liquid 
pressures: our experiment was at a pressure of 29 bar, while Korhonen's experi- 
ment was at 10.3 bar. Taking this pressure difference into account, the agreement 
is rather good. 

Generally, the height of the cf peak is not an intrinsic property of 3 He-B, 
since the cf-peak height depends in a nonlinear fashion on the field homogeneity 
and experimental geometry. In our measurements, qualitatively speaking, the cf- 
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Fig. 6 Thermal hysteresis of the critical texture transition fi c 2 between the parted 
and extended flare out textures: the plot shows the cf-peak at f! = 0.8rad/s with a vortex 
cluster {2 V = 0.1 rad/s, during a temperature sweep from high to low temperatures and back. In 
this cycle, the NMR line shape was recorded and the cf-peak was extracted after each spectrum 
was normalized in the reduced frequency domain. At high temperatures the order parameter 
was in the extended flare out texture and slowly cooled down. When the temperature cooled 
to T = 0.33 T c , the texture changed to the parted flare out texture. After reaching the low 
temperature point of T = 0.27 T c , the temperature sweep direction was reversed. After reaching 
a temperature of T = 0.354 T c as measured by the fork, the parted flare out texture flipped 
back to the extended flare out texture. Note that the relative noise of the cf-peak measurement 
increases in this figure with decreasing temperature due to the normalization of the spectra 
when the susceptibility decreases with decreasing temperature. 



peak height increases almost linearly in the parted flare out texture with increasing 
rotation velocity starting from fl c \ . At the texture transition from the parted to the 
extended flare out texture, the cf-peak height drops significantly when compared to 
the cf-peak height at the same rotation velocity in the parted flare out texture. The 
increase of the cf-peak in the extended flare out texture continues as a function of 
fi until its height slowly saturates, i.e. a maximum in absorption has been reached 
where all the area of the NMR line shape has moved "under" the cf-peak. As for 
the frequency shift Av of the cf-peak: in the parted flare out texture, the positive 
frequency shift from the Larmor frequency increases with increasing cf velocity. 
The frequency shift continues to increase in the extended flare out texture until it 
quickly saturates at Av = sin 2 j3 = | . Examples of the absorption and frequency 
shift measurements of the cf-peak have been published earlier and we refer to Fig. 
1 in [6] and Fig. 3 in [T6] , 

Kopu et al. pQ performed an extensive study on the cf-peak behavior and its 
dependence on the vortex number at high temperatures (T > 0.75T C ). Since the 
cf-peak height depends strongly on the exponentially decreasing density anisotropy 
in the low temperature limit T — > 0, in combination with the increasing textural 
transition velocity fi c i, the usefulness of this NMR measurement technique can be 
put into question. However, NMR measurements are still a useful technique down 
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Fig. 7 Critical textural transitions: figure shows the critical velocity i? c i (between simple 
and parted flare out textures) and J7 C 2 (between parted and extended flare out textures) given 
as the cf velocity Q c f = Q — Q v versus reduced temperature. The first critical velocity i? c i is 
constant down to T = 0.4T C and then increases rapidly. The error bars indicate the magnitude 
of the velocity step AQ when the transition was measured. The second critical velocity fl c i 
was measured in setup 2005 at constant rotation velocity by sweeping the temperature (T f: 
increase in temperature), and in the setup 2009 at constant temperature by varying rotation 
velocity (Q J [: rotation is increased by a step AQ). Error bars indicate the velocity step 
where the transition occurred. The pressure of the 3 He liquid was p = 29 bar. Our results are 
compared with data from Korhonen et al. (1990) at pressure p = 10.3 bar |24| . Larmor field 
values are given in the legend. 



to 0.2T C , when enough flow is applied (J7 c f > J? c i = 0.9rad/s). As an illustration, 
Fig. [8] shows a NMR line shape measurement of the vortex free state of 3 He-B at 
0.23 T c rotating at Q = 0.9rad/s: the cf-peak is significant. Calculations of spectra 
with varying sizes of vortex clusters are compared to this measurement in the plot 
and demonstrate that the cf-peak is still a sensitive measure of the number of 
vortices. 



8 Comparison of Measurements with Calculated Textures 

The measured NMR lines shapes converted in the frequency domain were fitted 
to numerical calculations using the same geometrical, environmental and NMR 
parameters as applicable for the experimental setup (see Section[5]). Measurements 
carried out during a continuous temperature sweep were accepted only if the rate 
of change of the temperature was within the absolute value \T\ < 5 • 10 -5 T c /s 
(includes stable temperature measurements) . The measured and calculated spectra 
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Fig. 8 Spectra as a function of equilibrium vortex cluster Q v ; The left figure shows 
the experimental (blue) and the optimal calculated (dashed green) line shape of vortex free 
supcrfluid 3 He at 0.232 T c rotating at Q = 0.9rad/s. The experimental line shape is an average 
of 3 measured spectra. The NMR resonance frequency was at u r { = 1.97 MHz. The optimal 
value for the field velocity parameter was Ajjv = 0.94kg/m 3 T 2 . Gray lines show the calculated 
line shapes with an increasing amount of vortices; we increased the equilibrium vortex cluster 
fi v with steps of AQ V = 0.05rad/s. The top right figure shows the frequency shift of the 
cf-peak as a function of the number of vortices in the central cluster, as expressed in terms of 
i? v and N . The bottom right figure show the extracted cf-peak height of the calculated NMR 
line shapes as a function the fl v and N. 



were plotted in the same figure on top of each other as seen in Fig.[T] We considered 
the fit to be good when (a) the NMR spectrum edge at high frequencies, which 
is determined by the B-phase longitudinal resonance frequency Hb(T, p), overlap 
with the edge in the measured spectrum; (b) the cf-peaks are nearly identical in 
height and frequency shift; (c) the difference in the normalized integrated area A A 
(see (|26[)) between the spectra is smaller than O.f (dimensionless in the reduced 
frequency domain). 

The fitting parameters are the temperature T, the field inhomogeneity AH/ H, 
and the field velocity parameter Ahv- The initial value for the temperature was 
taken from the fork reading, the field inhomogeneity from the NMR response in 
the Fermi-liquid state and the field velocity parameter Ahv with A(T) = A sc 
was calculated from (|9]) [7] . The fit using these initial values was typically rather 
poor: below T m 0.55T C the calculated cf-peak height was exceeding the measured 
cf-peak height and the measured frequency shift was less than expected. 

In the calculated spectra the frequency shift Av from the Larmor value is 
controlled by the temperature-dependent Leggett frequency Qb{T) and the cf- 
peak absorption by the field velocity parameter Ahv- The spread in the spectrum 
stems from the field inhomogeneity AH/H, which influences directly the cf-peak 
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Fig. 9 Longitudinal resonance frequency: Q-q(T,p) obtained from the fitted NMR line 
shape in the parted and extended flare out textures and a comparison with data from Ahonen 
et al. |28[ . The main panel shows the result of 321 runs (gray markers) and in binned form (blue 
circles). The inset shows the frequency shift from the Larmor frequency: in the parted flare out 
texture the value of the magnetic field was obtained from the intersection of a tangent of the 
edge and the base-line; in the extended flare out texture the magnetic field position of the 90° 
peak was used. The frequency of the cw-NMR resonance circuit was 965.0 kHz and the fluid 
pressure p = 29 bar. Fits were done using the formula (Q b /2tt) 2 = (1 -f 4 )(a-bf 4 +cf 6 ) ■ 10 10 
(reduced temperature T = T/T c , f2^/2n in Hz); with parameter values a = 8.73312, b = 
13.32121 and c = 1.51919 for our data; and parameter values a = 9.00301, b = 19.927 and 
c = 15.3442 for the data of Ahonen et al. (1976). 



height. We adjusted the temperature, which solved in first approximation the 
deviation in the frequency shifts from the Larmor value. By manually adjusting 
the field inhomogeneity and/or the field velocity parameter Ahv in the individual 
spectra did not give a unique solution for the best fit, since changes in these 
parameters are not orthogonal. As an illustration, see the cf-peak behavior for 
different values of Ahv in Fig. [TO] and consider that changes in AH/H also affect 
the peak broadening and peak height in linear NMR: at high counterflow velocities 
the cf-peak height increases with increasing Ahv as well as for decreasing values 
of AH/H. However, using both the parted and extended flare out line shapes at 
the same temperature when the texture transition happens at Q C 2 one can find a 
unique solution. Another method is to use the value of the critical velocity Q c i 
for the determination of Ahv: higher values of Ahv move il c i to a lower value. 
Subsequently a spectrum measured at high rotation velocity Q is used for the 
determination of the field inhomogeneity AH/H. After the use of either method, 
a next iteration by correcting the temperature improves the fit further. 
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Fig. 10 Uncertainty in the fitting procedure for the field velocity parameter Ajjv 
and isolation procedure of the surface peak: the main panel shows the measured spec- 
trum (blue) and the calculated spectrum (dashed green) using the optimal value for Ajjv- The 
different line shapes drawn with gray curves are calculated by varying Ajjv around the optimal 
value (the bottom inset shows the corresponding n- vector textures), from which the cf-peak 
dependence on Ajjv i s determined, which is used for estimating the uncertainty in Ahv> see 
1271 1 and (I28D . Taking the absolute difference between the measured spectrum (blue) and the 
calculated spectrum (dashed green) isolates the surface peak (red). Note that the area under 
the surface peak is the main contributor to the uncertainty in Xhv as it directly influences 
the normalization of the measured spectrum. 



The results of fitting the Leggett longitudinal frequency I?b as a function of 
the measured fork temperature are plotted in Fig. [9] and compared to data from 
Ahonen et al. [28]. The temperature reading for our data is taken from the fork 
resonance width while that in [28] was obtained by NMR measurement of fine 
platinum powder immersed in 3 He-B. The inset shows the frequency shift from 
Larmor, Av = v — v^, for both the parted and extended flare out textures: this 
result depends on the magnitude of the applied field H or rf frequency, while the 
result for the main panel is general. 

To obtain J?b directly from the measured frequency shifts using equations 
(|13 p ... (fl5|) . one has to resort to a feature with known angle /3. In the simple flare 
out texture, this would be the edge with the highest shift from the Larmor value; 
in the extended flare out texture the position of the 90°-peak can be used and 
in textures with high cf velocity the cf-peak position is at AD = 0.8. See for 
example Fig. 1 in [24_. Each method has its difficulty in determining the correct 
value, in particular when the position of the observed feature is on a slope. From 
numerical calculations we see that the observed 90°-peak is shifted due to the field 
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inhomogeneity. The error in the frequency shift is estimated to be below 2%, which 
is considered large. The error can be decreased by measuring spectra with a lower 
value of ll>l such that the spectrum is wider, see (|13p . but this also reduces the cf 
peak sensitivity. 

Even for the best fit there still exists a small discrepancy between the measured 
(Xm) and calculated (Xc) line shapes, which is expressed by the absolute value of 
the difference between the two line shapes normalized to the area of the measured 
spectra, 

AA= S\x^)-M\d r 

We use this value for the uncertainty in Ahv by calculating the spectra with 
different values of Ahv around the optimum, which gives the relationship between 
the cf-peak height Xcfp an d Ahv, i.e. Anv(Xcfp)- By defining the uncertainty in 
the cf-peak height Xcfp±> the uncertainty in Ahv is also found 

Xcfpi = Xcfp • (1 ± AA), (27) 
A H v± = A H v(Xrf P ±)- (28) 

The procedure is demonstrated in Fig. 1101 In that figure, the optimal value of 
Ahv for fitting the calculation and the measurement is Ahv = 5.3kg/m 3 T 2 , which 
gives a cf-peak height x'cfp = 0.053 kHz -1 . Increasing the field velocity parameter 
to Ahv = 6.3kg/m 3 T 2 gives the first gray line above the optimum calculated 
spectrum and the cf-peak increases to x'cfp = 0.063 kHz -1 . At the upper bound 
value of Ahv = 5.8kg/m 3 T 2 , the absolute value of the difference between the 
line shapes, A A, is fully attributed to the uncertainty in Ahv- In this figure the 
relation between Ahv and x'dp is approximately linear, but this is not generally 
the case. In particular, when the total area under of the NMR line shape is almost 
completely under the cf-peak (high counterflow) , the dependence of Xcfp on Ahv 
becomes almost constant and this increases the uncertainty Ahv±- 

Fig. [TT] shows the measured field velocity parameter Ahv as a function of 
inverse temperature as measured with the forks. Gray markers are the results of 
individual measurements and blue markers are averages of this data within a bin 
of inverse temperature width A(T C /T) = 0.25. In the high temperature regime, 
our data compares well with those of Korhonen et al. |25j . In the low temperature 
regime additionally data is included from our measurements of Ahv using coherent 
precession NMR mode. This mode corresponds to the Bose-Einstein condensation 
of magnons on the ground level of a 3-dimensional potential well, which is created 
by the order parameter texture and applied profile of the static NMR field [26] . 
The position of the level, which can be measured with NMR, is a sensitive probe 
for the dependence /3(r) close to the axis of the sample (while the linear NMR 
response, explained above, is more sensitive to the f3{r) behavior at larger radii). 
The measured slope 4^ | r= o can be connected with the relevant textural parameters 
via numerical calculations of the texture. The technique is explained in |27j . where 
measurements are done in the equilibrium vortex state and the vortex textural 
parameter A is extracted. In this work we have performed analogous measurements 
in the vortex-free state using the upper spectrometer in the 2005 setup and fitted 
the magnon level positions using Ahv as a fitting parameter. The results obtained 
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Fig. 11 Field velocity parameter Ahv as a function of inverse temperature. Our 

data using NMR line shape measurements (gray markers, 350 individual NMR spectra), Ahv 
was obtained using the fitting procedure as described in Fig. 1101 The error-bars indicate the 
boundaries of Ahv as determined by the area mismatch, see H28I I. Blue markers show averages 
within the inverse-temperature bins; on these markers the error-bars indicate ±1 standard 
deviation. In the low temperature regime, data is shown where Ahv is obtained from fitting 
Ahv m experiments with magnons in the potential well. In the high temperature regime data 
is shown from Korhonen et al. 25 . Solid and dashed lines represent the model of the field 
velocity parameter Ahv ® with the energy gaps A BC 1241 ) and Zi sca i ct i (1251 . 



using this non-linear NMR mode agree well with the measurements using linear 
NMR response. 

Calculated values for Ahv (see (|9])) are shown in Fig. [TT] for the energy gap 
with strong-coupling correction, A SC (T) (dashed line) [7], sometimes also referred 
to as weak-coupling plus (A wc +), and for the scaled energy gap A sca \ e d(T). For 
the definitions of these energy gaps, see (|24[l and (|25|) . Measurements at pressures 
p < 19.4 bar by Davis et al. using transverse and longitudinal sound in 3 He-B 
suggest the value of the energy gap A sc . 

The temperature in our data is obtained from the resonance frequency width 
Au of the fork by sweeping the fork around its resonance frequency. The resonance 
width at high temperatures (T > 0.34T) is then calibrated against the MCT as 
a function of temperature. In the low temperature regime the dependence on 
temperature is extrapolated as an exponential dependence on the energy gap: 
Auj oc exp(— A/T). All temperature measurements of our data in this paper use the 
energy gap value zA ex p(0) = 1.97 fee T c . When we compare the experimental data 
with predictions of the theory using alternative energy gaps (|23|) . . . (f25j) . we have to 
rescale the temperature measurement of the experiments, i.e. use the energy gap 
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Fig. 12 Surface peak measurements: the main panel shows the reduced frequency shift 
Aii = sin 2 /3 of the 90°-peak (squares) and the surface peak (stars) as a function of reduced 
temperature. Gray markers are individual measurements while blue markers are averages in 
temperature bins of width AT = 0.015 T c . The inset shows the reduced textural energy U(r) = 
sin 2 /3(r) as a function of reduced radius at different temperatures in a cylinder of radius 
R = 3 mm. The curves illustrate the potential well at the cylindrical surface, which is formed 
between the location from where the 90° -peak originates and the surface boundary at r = R. 



A consistently. When we rescale the temperature of our data in Fig. [TT] using the 
strong-coupling corrected gap A sc , data points above T c /T « 3 move to higher 
inverse temperature with maximum shift at T c /T = 5. However, the temperature 
shift is minimal and the measured Ahv has still a better correspondence using the 
scaled energy gap with the low temperature limiting value Zi e xp(0) = 1.97/cbTc. 



9 Surface Spin Resonance 

In our measurements of the extended flare out texture at temperatures T < 0.45 T c , 
three absorption peaks have been observed in the NMR spectra (see Fig. [JJ) . The 
origin of the cf-peak and the 90°-peak is identified in terms of the calculated 
texture and the local NMR response. The cf-peak results from a large region of the 
sample where the angle ft is (approximately) constant and the 90°-peak from the 
maximum reduced frequency shift when (3 crosses 90°, i.e. AD = sin 2 (90°) = 1. In 
Fig. [10] the absolute value of the difference (red line) is taken from the measured 
and calculated spectra which isolates a third peak. We call it the surface peak 
based on the interpretation described below, which is an extension of one of the 
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brief suggestions by Salomaa |14| . A depiction of the peak appears also in later 
publications, see Fig. 1 in [23] and Fig. 4 in [31] . 

Since the new peak is only observed in the extended flare out texture, we 
suggest that it is associated with the particular n(r) distribution close to the 
cylinder wall. In Fig. [12] the reduced frequency shift Av of the surface peak is 
plotted against the reduced temperature T/T c . The result shows that the shift 
from the Larmor value decreases with decreasing temperature. The figure also 
contains the reduced frequency shift of the 90°-peak. When the separation of the 
surface peak from the 90° peak is small, it becomes hard to distinguish these 
two peaks due to their overlap. The inset of Fig. [12] shows the reduced textural 
energy U(r) = sin 2 /3(r) as a function of reduced radius at different temperatures. 
The decreasing energy at large radii creates a potential well for magnons between 
/3(r) = 90° (AD = sin 2 f3 = 1) and the cylindrical surface at r = R. We suggest 
that this potential well is the origin of the surface peak. The magnon eigenstates 
in the potential give rise to a coherent non-local spin- wave resonance [32]. We 
suggest that the absorption at the frequency of the ground state forms the surface 
peak. With decreasing temperature the texture produces a potential well with 
lower energy, which corresponds qualitatively to a smaller frequency shift from 
the Larmor frequency. 

10 Conclusions 

In this work we have compared measured NMR line shapes of rotating super- 
fluid 3 He-B in a cylindrical environment with calculated line shapes using the 
Hydrostatic Theory of 3 He-B. Quantitative agreement can be achieved by ad- 
justing the field velocity parameter Ahv- The measured values of the parameter 
Ahv show a deviation from the prediction for A(T) with strong-coupling correc- 
tion, i.e. Zi sc (0) = 1.87fcBT c . A better correspondence between the measured and 
predicted values of Ahv is established when A(T) is shifted toward a higher value 
^exp(O) = 1.97 k&Tc as earlier measured by Todoshchenko et al. (2002). 

Measurements on the susceptibility xb/xn have been extended to lower tem- 
peratures (T = 0.2 T c ). The low temperature limiting value of xb/xn agrees 
with the currently accepted value of the anti-symmetric Fermi-liquid parameter 
Fq = —0.75 at 29 bar liquid pressure. 

By fitting the calculated and measured NMR spectra, we can determine the 
B-phase longitudinal resonance frequency J?b (T, p) , which shows a slight deviation 
below T < 0.5 T c from the data measured by Ahonen et al. [28]. Due to uncertainty 
in the measured pressure, we can not claim our result in the pressure dependent 

(T, p) is better. Since we here fit the whole spectrum which includes the broad- 
ening due to the field inhomogeneity, we believe that this procedure represents an 
improvement over the method used by Ahonen et al. 

Taking the absolute difference between the measured and calculated spectra re- 
vealed in the extended flare out texture below 0.45 T c a new absorption peak, which 
we interpret as a spin-wave resonance in the potential well between the 90°-peak 
and the cylindrical surface. There is qualitative agreement between the reduced 
frequency shift of this peak and the shift in the potential well with decreasing 
temperature. 



21 



R. de Graaf et al. 



Acknowledgements The authors thank Juha Kopu for making his optimization program 
available. Erkki Thuneberg and Grigory Volovik are thanked for the theoretical discussions. 
Contributions to the measurements and the experimental setup by Antti Finne, Roman Sol- 
ntsev and Rob Blaauwgeers are greatly appreciated. This work is supported in part by the 
Academy of Finland (Centers of Excellence Programme 2006-2011) and EU 7th Framework 
Programme (grant 228464 Microkclvin) 



References 

1. J. Kopu, R. Schanen, R. Blaauwgeers, V.B. Eltsov, M. Krusius, J.J. Ruohio and E.V. 
Thuneberg, J. Low Temp. Phys. 120, 213 (2000) 

2. R. Hanninen, V.B. Eltsov, A. P. Finne, R. de Graaf, J. Kopu, M. Krusius and R.E. Solntsev, 
J. Low Temp. Phys. 155, 98 (2009) 

3. R.E. Solntsev, R. de Graaf, V.B. Eltsov, R. Hanninen and M. Krusius, J. Low Temp. Phys. 
148, 311 (2007) 

4. R. de Graaf, R. Hanninen, T.V. Chagovets, V.B. Eltsov, M. Krusius and R.E. Solntsev, J. 
Low Temp. Phys. 153, 197 (2008) 

5. V.B. Eltsov, A.I. Golov, R. de Graaf, R. Hanninen, M. Krusius, V.S. L'vov, R.E. Solntsev, 
Phys. Rev. Lett. 99, 265301 (2007) 

6. V.B. Eltsov, R. de Graaf, P.J. Heikkinen, J.J. Hosio, R. Hanninen, M. Krusius, and V.S. 
L'vov, Phys. Rev. Lett. 105, 125301 (2010) 

7. E.V. Thuneberg, J. Low Temp. Phys. 122, 657 (2001) 

8. J. Kopu, J. Low Temp. Phys. 146,47 (2007) 

9. D.D. Osheroff and W.F. Brinkman, Phys. Rev. Lett. 32, 584 (1974) 

10. R. Balian and N.R. Werthamer, Phys. Rev. 131, 1553 (1963) 

11. D. Vollhardt and P. Wolfle, The Superfluid Phases of Helium 3, 81, Taylor & Francis, 
L ondon (1990) 

12. http://www.netlib.org/opt/tn 

13. P.J. Hakonen, M. Krusius, M.M. Saloma, R.H. Salmelin, J.T Simola, A.D. Gongadze, G.E. 
Vachnadze and G.A. Kharadze, J. Low Temp. Phys. 76, 225 (1989) 

14. M.M. Salomaa, J. Phys-Condens Matter 2, 1325 (1990) 

15. W.F. Brinkman, H. Smith, D.D. Osheroff and E.I. Blount, Phys. Rev. Lett. 33, 624 (1974) 

16. V.B. Eltsov, R. de Graaf, P.J. Heikkinen, J.J. Hosio, R. Hanninen and M. Krusius, J. Low 
Temp. Phys. 161, 474 (2010) 

17. D.S. Greywall, Phys. Rev. B 33, 7520 (1986) 

18. T.D.C. Bevan, A.J. Manninen, J.B. Cook, H. Alles, J.R. Hook and H.E. Hall, J. Low 
Temp. Phys. 109, 423 (1997) 

19. LA. Todoshchenko, H. Alles, A. Babkin, A.Y. Parshin and V. Tsepelin, J. Low Temp. 
Phys. 126, 1449 (2002) 

20. E.R. Dobbs, Helium Three, 512, Oxford University Press, New York (2000) 

21. H.M. Bozler, P.D. Saundry, Inseob Hanh, S.T.P. Boyd and CM. Gould, J. Low Temp. 
Phys. 89, 5 (1992) 

22. H.N. Scholz, Ph.D. thesis, Ohio State University, 1981 (unpublished) 

23. Inseob Hahn, S.T.P. Boyd, H.M. Bozler and CM. Gould, Phys. Rev Lett. 81, 618, 1998 

24. J.S. Korhoncn, A.D. Gongadze, Z. Jami, Y. Kondo, M. Krusius and E.V. Thuneberg, 
Phys. Rev. Lett. 65, 1211 (1990) 

25. J.S. Korhonen, Y.M. Bunkov, V.V. Dmitriev, Y. Kondo, M. Krusius, Y.M. Mukharskiy, 
U. Parts and E.V. Thuneberg, Phys. Rev. B 46, 13983 (1992) 

26. Yu.M. Bunkov and G.E. Volovik, Phys. Rev. Lett. 98, 265302 (2007) 

27. V.B. Eltsov, R. de Graaf, M. Krusius and D.E. Zmeev, J. Low Temp. Phys. 161 (2011) 
[DOI: 10.1007/sl0909-010-0285-l]. 

28. A.I. Ahonen, M. Krusius and M.A. Paalanen, J. Low Temp. Phys. 25, 421 (1976) 

29. D.D. Osheroff, W.J. Gully, R.C Richardson and D.M. Lee, Phys. Rev. Lett. 29, 920 (1972) 

30. J. P. Davis, H. Choi, J. Pollanen and W.P. Halperin, J. Low Temp. Phys. 153, 1 (2008) 

31. A.P. Finne, S. Boldarev, V.B. Eltsov and M. Krusius, J. Low Temp. Phys. 136, 249 (2004) 

32. Yu.M. Bunkov, V.B. El tsov, R. de Graaf, P.J. Heikkinen, J.J. Hosio, M. Krusius and G.E. 
Volovik. [arXivTl 002.1674vl [cond-mat.quant-gas], (2010) 



